Phase transitions and ordering of confined dipolar fluids 
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We apply a modified mean-field density functional theory to determine the phase behavior of 
Stockmayer fluids in slitlike pores formed by two walls with identical substrate potentials. Based on 
the Carnahan-Starling equation of state, a fundamental-measure theory is employed to incorporate 
the effects of short-ranged hard sphere - like correlations while the long-ranged contributions to 
the fluid interaction potential are treated perturbatively. The liquid-vapor, ferromagnetic liquid - 
vapor, and ferromagnetic liquid - isotropic liquid first-order phase separations are investigated. The 
local orientational structure of the anisotropic and inhomogeneous ferromagnetic liquid phase is also 
studied. We discuss how the phase diagrams are shifted and distorted upon varying the pore width. 

PACS numbers: 



I. INTRODUCTION 

The structure and thermodynamics of confined dipo- 
lar fluids (such as molecular liquids, ferrofluids, and elec- 
trorheological fluids) have recently attracted considerable 
attention. Theoreticali'^'^, computer simulation^i^i^iiiSii^, 
and experimentaliSiii studies have provided significant 
information about the orientational and spatial arrange- 
ment of dipolar particles in the vicinity of solid walls. The 
strongly anisotropic and long-ranged character of dipolar 
forces causes particular difficulties for the theoretical de- 
scription of these systems. As long as one is not aiming 
for a quantitative description of a specific system but for 
general phenomena and trends, the so-called Stockmayer 
model has turned out to be rather usefuli^. It considers 
spherical particles interacting with Lennard- Jones (LJ) 
pair potentials plus pointlike permanent dipoles at their 
centre. In this paper we adopt the magnetic language, 
assuming that the particles carry magnetic dipole mo- 
ments of strength m because the main applications we 
have in mind are ferrofluids. (The results are identical 
for ferroelectric particles.) 

Bulk dipolar fluids such as the Stockmayer model sys- 
tem may exhibit three distinct fluid phases: isotropic 
vapor, isotropic liquid, and ferromagnetic liquidi^i^. (In 
the present study we focus on sufficiently high temper- 
atures so that freezing is not a concernii.) Accord- 
ingly the following first-order phase transitions can oc- 
cur: isotropic vapor - isotropic liquid, isotropic liquid - 
ferromagnetic liquid, and isotropic vapor - ferromagnetic 
liquid. The phase diagrams exhibit the usual liquid-gas 
critical point, a triple point, a tricritical point, and a 
line of critical points corresponding to second-order phase 
transitions between the isotropic and the ferromagnetic 
liquid. For bulk dipolar fiuids these phase transitions 
and critical points and lines have been detected in vari- 
ous theoretical studie a-^^i-^^'-^'^i-^^'-^^ . 

However, to the best of our knowledge there is no theo- 
retical study which predicts the complete global phase di- 
agrams and the structure of confined dipolar fiuid phases. 



Recently Gramzow and Klapp^ have studied the phase 
behavior and orientational structure of Stockmayer fluids 
confined to slit pores. Within their modified mean-field 
(MMF) density functional theory (DPT) they considered 
only the subspace of spatially homogeneous local densi- 
ties throughout the slit pore and focused on the orien- 
tational structure only. This allowed them to study the 
infiuence of the confinement on the phase behavior but 
neglecting the inevitable spatial inhomogeneities of the 
phases under study. Our results will show that in the 
case of confined dipolar liquids, as for other fiuids, there 
are strong spatial inhomogeneities in the vicinity of the 
walls. This inhomogeneous character of the liquid phases 
becomes more pronounced upon increasing of the density. 
These phenomena require a more sophisticated theory in 
order to describe the structure and the phase equilibria 
of confined dipolar systems. 

In the present study we also use the modified mean- 
field (MMF) DPT approximation. However, for the de- 
scription of the hard sphere reference system, instead 
of the homogeneous bulk Carnahan-Starling free energy 
expression, a fundamental measure theory free energy 
functionalii is applied in order to capture the short- 
ranged hard sphere - like correlations. Our results for 
the confinement induced shifts and the distortions of the 
phase boundaries relative to the bulk ones are compared 
with the corresponding Kelvin equation. 

The paper is organized as follows. In Sects. 2 and 3 
the model and the extension of the MMF density func- 
tional theory for confined dipolar fluids are presented and 
developed, respectively. In Sec. 4 the calculation of the 
orientational order and the physical meaning of the or- 
der parameters are discussed. Sections 5 and 6 present 
the details of the Euler-Lagrange and of the Kelvin equa- 
tions, respectively. The results of our calculations and a 
discussion are given in Sec. 7. Certain important com- 
putational details are summarized in the Appendices A, 
B, and C. 
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II. MICROSCOPIC MODEL 



We study so-called hard core Stockmayer fluids which are characterized by the interaction potential 



Us{vi2,UJi,UJ2) 



00 



ri2 < cr, 



(1) 



r 



where 



= 3[mi(wi) • ri2][m2(cJ2) -^12] - [mi(wi) • 1112(^2)] (2) 

is a rotationally invariant function. In Eq. ([T]) the first 
term is the Lennard- Jones potential with length and en- 
ergy parameters a and e, respectively. The second term 
in Eq. ^ together with Eq. ^ describes the dipole- 
dipole interaction potential, where particle 1 (2) with 
diameter a is located at ri (r2) and carries at its center 
a point dipole moment of strength m with an orienta- 
tion given by the unit vector mi(a;i) (m2(w2)) with po- 
lar angles wi = (6*1, <^i) {uj2 = (6*2,02)); ri2 = ri - r2 
is the difference vector between the centres of particle 
1 and 2, ri2 = |ri2|, and fi2 = ri2/ri2 is a unit vec- 
tor with orientation W12 = {0i2,4>i2)- In order to apply 
density functional theory this interaction potential is de- 
composed into a short-ranged repulsive part 



walls we choose the external substrate potential 



Mr(fl2) = Uhsiri2) 



, ri2 > cr 
00 , ri2 < cr 



and into a long-ranged excess part 

We2-c(ri2,t^i, W2) = Q{ri2 - a)us{ri2,uji,uj2). 
Thus the total pair potential is given by 

Ms(ri2, wi, W2) = Uriri2) + Ue:cc(ri2,wi,a;2). 



(3) 



(4) 



(5) 



This decomposition lends itself to choose a hard-sphere 
reference system characterized by u, (fi2) = Uhs{r 12) for 
which reliable approximations for the free energy are 
known. The excess part Ue^c(ri2, wi, 0)2) will be treated 
perturbatively in an appropriate way. 

For the fluid-wall interaction we consider two different 
potentials. In the case of a purely hard wall we use the 
hard-wall potential 



Uwh{z) 



, Izl < {L-a)/2 
00 , |z| > (L-cr)/2, 



(6) 



where z is the coordinate normal to the walls and L is the 
distance between the surfaces of the hard, parallel walls, 
which we call repulsive walls. For the study of attractive 



2tt , 



■L/2+, 



^wa (^) 



, \z\<{L-a)/2 
, \z\>{L-a)/2, 
(7) 



where sets the energy scale. This equation can be ob- 
tained from the assumption that both walls consist of LJ 
particles interacting with the fluid particles of the same 
size also via the LJ potential with the energy scale e^. 
Averaging the LJ interactions over all positions of the 
wall particles and approximating the averaged repulsive 
part of the LJ potential by a hard core potential one ob- 
tains the fluid particle - wall interaction potential given 
by Eq. ^. 



III. MODIFIED MEAN-FIELD DENSITY 
FUNCTIONAL THEORY FOR 
ONE-COMPONENT FLUIDS 



For a one-component nonuniform fluid configuration 
p(r, uj) denotes the number density of dipolar particles 
at a point r ~ [x, y, z) with an orientation uj = {6, (j)) 
relative to a spatially fixed coordinate system. The total 
number density of particles, independent of orientation, 
is given as 



/f2iT fir 
dLop{r,uj) ~ d(j) d9 sin{9)p{r,iu). 
Jo Jo 



(8) 



This allows one to split /5(r, uj) into the total number den- 
sity /5(r) and a normalized space- and angle-dependent 
orientational distribution function a{r,uj): 



p{r,uj) = p(r)a(r,w). 



1. 



(9) 



Within density functional theoryi^ the equilibrium den- 
sity distribution p(r,u!;T, ^) of an inhomogeneous fluid 
in the presence of an external potential minimizes the 



3 



grand canonical potential 

n[{p{v,Lo)},T,^i] 

^F[{p{v,io)],T]+ [ d^rdLop{Y,io){u^{r,uj)^ p) 



+ (rrK{v)\l- duja{v,uj)], (10) 



and thus satisfies the Euler-Lagrange equation 



5n[{p{v,w)},T,p] 



5p{r,u)) 



(11) 



p[v .uj)— p[v .UJ\T 



Here F is the Helmholtz free energy density functional of 
the system and p is the chemical potential. The last term 
in Eq. (|10p takes into account the normalization condi- 
tion (Eq. ^) for the orientational distribution function 
(5;(r,w) by a Lagrange parameter k.{v). Due to the two 
distinct contributions to the particle interaction poten- 
tial (Eq. (O) the free energy functional F decomposes 
into an ideal gas term F"^ and two corresponding parts: 



F[{p{v,u)},T]^F^%p{v,u:)},T] 
+F-^[{/5(r)},T]+^^-[{p(r,c.)},r], 



(12) 



where F'^'^^ is the reference and F'^^'^ is the excess free 
energy density functional. 

In order to describe the planar wall - fluid interfaces, 
we consider a slab-shaped macroscopic system with the 
surfaces of the slab parallel to the x — y plane. The 
distance between the parallel surfaces is L. Thus, apart 
from spontaneous symmetry breaking like freezing, the 
equilibrium configuration of the system is inhomogcneous 
only in the z direction and translationally invariant in the 
X and y directions. Under these circumstances the total 
number density p{v) = p{z) is a function of z only and 
the orientational profile a(r,w) = a{z,uj) depends on z 
and the polar angle lo. In a coordinate system fixed in 
space the actual orientational distribution function 



following m = —m.) The expansion coefficients aim{z), 
which can be interpreted as orientational order parame- 
ters, are related to the full distribution function via 



aimiz) = / duja{z,u)Yi*^{u). 



(14) 



As presented below, for a slab-shaped system of transver- 
sal size L one obtains explicit expressions for the three 
terms of the Helmholtz free energy density functional (see 
Eq. (HH)). 

A. The free energy functional for the ideal gas 

For the slablike shape of the sample under considera- 
tion the ideal gas contribution has the form 

F"^^^ dz/5(z)[ln(p(z)A3)-l] 



dzp{z) I (iijjd(z, a;) ln[47rd(z, w)], (15) 



where A is the lateral cross-sectional area of the system, 
A is the de Broglie wavelength, and (3 = l/{kBT) is the 
inverse temperature. 



B. The reference free energy functional 

For the description of the hard sphere reference free en- 
ergy functional we adopt the fundamental measure the- 
ory, which was initially proposed by Rosenfel d^^'^° : 



A /-^/^ 



PJ- 



L/2 



dz^{{na{z)}); 



(16) 



the function $ and the weighted densities na{z) are de- 
fined in Appendix A. 



a{z,Lj) 



E 



almiz)Ylr,^{uJ), 



(13) 



can be expanded in terms of the spherical harmonics 
Yim{uj) where aoo = l/Vin due to the normalization 
and = (— l)™aim because a is real. (Here and in the 



C. The excess free energy functional 

We approximate the excess free energy functional 
in terms of the modified mean-field density functional 
theory2ii^i2^. Using Eq. (fT^ for the present slablike 
geometry one obtains the following expression: 



" = ■5~T^ / '^^i / / dxidyi I dx2dy2 x 

Stt^PJ_l/2 J -l/2 J a J a 

p{zi)p{z2) ^ ^ diimi(^l)di2m2(^2) X 

dw^ / dc^2>i„„,(a;i)e-"-^('--)(l-e-""=('--^"i^-^))yi,™,(w2). 



(17) 
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Here and in the following the summations over li, mi are 
taken as 



OO li 



(18) 



L=0 mi=-li 



In order to extract the explicit proportionality to the 
large cross-sectional area A we introduce the sum and the 
difference of the lateral coordinates and {x2,y2)- 



a:i2=xi-X2, 2/12 2/2, 

xi +X2 yi+ 2/2 

— — ' y = — — ■ 



(19) 



In terms of these new variables the integration of any 
function f{r 12) can be written as 



dxidyi / dx2dy2f{ri2) 
J A 



dxdy / dx i2dy 12 f{r 12) ^ A dxi2dyi2firi2){20) 
A J A J A 



In order to proceed we replace the rectangular coordi- 
nates (a;i2,?/i2) by polar coordinates {Ri2,4>i2) and take 
into account Eqs. ((T7]) and (PO)) for the ea;cess free energy 
functional: 

jpexc r^/'2 

—— = dzi dz2 p{zi)p(z2) X 

^ J-L/2 J-L/2 

(Z1-Z2), (21) 

where 

1 f°° 

Ulimihm.2{zi - Z2) = --5-^ / dRi2Ri20{ri2 - (t) X 



12 / dlVi I du)2Yl^miiL^l)fMiri2,UJl,U)2)Yi^m2i^2)- 

(22) 



In Eq. (|22p fi\i denotes the Mayer function of the excess 
potential: 

/M(ri2, 0.1, 0.2) = e-'3"=-('--''^-'^^) - 1, (23) 

where in accordance with the planar polar coordinates 
used for the slab-geometry 



ri2 ^ \ Rh + {zi - Z2y 



(24) 



Via the vector ri2 the rotationally invariant function 
D{uji,oj2,oji2) (Eq. ([2])) depends on 1912 as well. There- 



fore in Eq. ([221) the substitutions 

R12 



sinz9i2 = ■ 

C0Sl?i2 = 



Vi??2 + (^1-^2)2' 
Zl - Z2 



(25) 



have to be carried out prior to the integration over the 
variable -R12. In our calculation we have expanded the ex- 
ponential term in Eq. (|23p . containing the dipole-dipole 
interaction, into a Taylor series: 



-1 + e 



-PnL.j(.ri2) 



1-t- 



/3m2 ^ , 1 I3^m\^^ 



{ri2f 2!(ri2)' 



(26) 

In accordance with Eqs. and (|26p ui-^rnii2m2 is given 
as a Taylor series, too: 

_ (0) 

Ulimil2m2 — '^Iimil2m2 

+u\^^ , x(m2)i+up) , X (m2)2..., (27) 

Iimil2rn2 ^ ' Iimil2rn2 V / ' \ ' / 

with the coefficients 



xe(i?f2 + z -CT^)- 



*iimi/2m2 

exp 



f3 



i-l fOC 



2*! Jo 

pULJ (v/i??2 + 



di?12-Rl2 



{R\^ + z2)3»/2 

(0 ( ^ 



Iimil2m2 



In Eq. (f28|) the coefficients ^i'm^i^mj can be expressed 
as 



(28) 



(i(?!)i2 y dtJi j dijj2Yi^,m{i^i) 
X [i:»(tJi2, wi, a;2)]*y(2m2 (a;2)- 



(29) 



In Appendix B we present the nonzero coefficients 



A 



(i) 

Iimil2m2 



for i, li, I2 = 0, 1, 2. In the bulk limit the ex- 



pressions for ui^jn-^i^jn^ reduce to the corresponding ones 
in Ref J^: 



lim 

L — ^00 



L/2 



dz u 



-L/2 



(i) 

litnil2Tn[ 



iz) 



r, (30) 

where the coefficients for li = 0,1,2,3,4, are given 
by Eqs. (4.9)- (4. 13) in RefJ^. It is important to note 
that among the coefficients Ui-^rnii2m2 there are terms 
which arc nonzero in a slablike geometry, but which van- 
ish in the bulk limit. 
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IV. MAGNETIZATION AND ORDER 
PARAMETERS 



so that the modulus Adxy of the in-plane magnetization 
is 



Focusing on the structure of confined dipolar fluids, 
here we are particularly interested in the occurrence of a 
spontaneous magnetization. In the slab geometry, which 
can be considered as the limiting case of an oblate ellip- 
soidal sample with vanishing aspect ratio, the demagne- 
tization field vanishes if the fluid sample is magnetized 
along a spontaneously chosen direction within the xy 
plane. For such a type of magnetization the formation 
of various domains can be excluded. Therefore the con- 
fined fluid with the magnetization within the xy plane 
is comparable to a single domain bulk fluid in a needle- 
like volume with longitudinal magnetization. The free 
energies of these fluid samples can be mapped onto each 
other—. The definition of the local magnetization is 



(31) 



M(z) = m/9(z) / duja{z,uj)m.{uj). 



By inserting the expansion in terms of the spherical har- 
monics and performing the integration we obtain 



47r 



1/2 



-\/2i?e[aii(z)] 



M(z) - ( Y j T^p{z) <j +V2/m[«n(^)] ) ■ (32) 

aw{z) 

For Mz ~ 0, i.e., aio(z) = 0, the magnetization within 
the xy plane is measured by the order parameter 



1/2 



((i?e[«n(z)])2 + (/m[an(^)])')'^' 



3-; 



\oiii{z)\, 



(33) 
I 



(34) 



The higher coefficients of the oricntational distribution 
function a{z, lu) (Eq. (fT3|) ) may also be considered as lo- 
cal oricntational order parameters. Nonzero coefficients 
for I = 2 indicate a certain type of ordering of the dipoles. 
An especially interesting quantity is the order parameter 
a2o{z) that describes the orientation of dipole-axes rela- 
tive to the z axis. This is related to the zz component of 
the local quadrupole tensor as 



.(^) 



duja{z, uj)P2{cos 6) 



1/2 



020(2), 



(35) 

where P2 is the second order Legendre polynomial. 
Qzz ~ 1 corresponds to the perfect ordering of dipole 
axes normal to the walls while Q^z — ~l/2 indicates 
the perfect ordering along directions in the xy plane. 
Randomly oriented dipoles correspond to Qzz = 0. 



V. EULER-LAGRANGE EQUATION 



The equilibrium configuration characterized by p{z) 
and a(z, uj) for given values of T and p follows from min- 
imizing the total grand canonical functional with respect 
to p{z) and the function a{z,uj): 



(3 5Vt 



A5p{z) 



p—p^ a— Q 
L/2 



ln/)(z) + [{p(z)}] + / doj Q!(z, w) ln[47ra(z, w)] 



+2[i ai^,m{z) / dz' p{z')ui^rmi2mAz - z')ai^rnAz') ~ Pt^- + I3uw{z) = 

I ; J-L/2 



and 



/3 sn 



A Sa{z, Lo) 



= p{z) (1 + ln[47ra(z, uj)]) - I3k{z) 



p—p, a— a 



.L/2 

+2f3p{z) / dz'p{z') V ai,,m{z') V u/i^ij^m^ (z - z')y,*„, (w) = 

J-L/2 , — , — 



(36) 



(37) 



where t\J[{p{z)} 



5F"^[{p(z)}] . 



, . 1 r- ''° 1 /c ; J ^^(f) 1 i- r i- \ r corresponding set of integral equations is discussed in Ap- 

tial functional (nrst order direct correlation function) of ^ ° ° ^ ^ 



is the chemical poten- the hard sphere reference system (see Appendix A) . The 
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pendix C. As solutions of Eqs. p6|) and ([37]) p{z;T,^) 
and a{z, w; T, /i) become functions of T and /i. 



VI. PHASE EQUILIBRIA AND KELVIN 
EQUATION 

The phase coexistence curves fJ.(T) and the coexisting 
densities and orientational order parameters follow from 
requiring the equality of the grand potentials for the co- 
existing phases / and //: 



n[{p^'\z; T, p), a^'\z, u; T, p)}, T, p] = 
n[{p^''^ (z; T, p),a^''^ (z, u;; T, ^i)}, T, p] 



(38) 



The functions pW(z;T,/x) and a^'\z,uj;T, ^l) {i = I, II) 
denote the corresponding equilibrium density and orien- 
tational distribution functions obtained from Eqs. (|36p 
and (j37p . Equation (|38p requires simultaneous solutions 
of the Eulcr-Lagrange equations and (|57|) for a wide 
range of temperatures and chemical potentials. As men- 
tioned in Sec. 1, we consider three kinds of two-phase 
equilibria: isotropic liquid - isotropic vapor, isotropic 
vapor - ferromagnetic liquid, and isotropic liquid - fer- 
romagnetic liquid. In order to estimate in leading order 
of 1/L the location of the phase coexistence curves, the 
Kelvin equation (see, e.g. Refi^ and references therein) 
is used to obtain the chemical potential difference be- 
tween the bulk phase points and the corresponding ones 
for slabs of thickness L at a given temperature and for 
two identical confining walls: 



A^i = n[L) - p[co) = 

^ Pi - Pi I 



2 7/,// cos 9/ 
L pi 



Pii 
(39) 

where pi [i = I, II) arc the bulk densities of the coex- 
isting phases, 7^,^ arc the corresponding wall-fluid inter- 
facial tensions, 77,77 is the surface tension between the 
coexisting fluid phases / and //, and 87 = arccos[(7u,77 — 
7-!«7)/77.77] is the contact angle of a drop of fluid / with 
the wall (Young's equation). In the slab-shaped system, 
for a fluid which is in contact with two identical walls, the 
surface tension ^^ji follows from the equilibrium density 
p{z;T,p) and oriental distribution a{z,u);T, p) profiles 
as 

Iwi = 7; 1™ {Lp + flL[{p{z\T,^i),a{z,uj;T,p)},T,p]), 

(40) 

where p is the bulk pressure and fi^ is the grand poten- 
tial per cross-sectional area A of the slablike system of 
transversal size L. 



VII. RESULTS AND DISCUSSION 

In the following we shall use dimensionless quanti- 
ties: T* = ksT/t as reduced temperature, p* = (p — 
ksT \n{A^ /a^))/e as reduced chemical potential, p* ~ 



77 /-L/2 dzp{z) as reduced averaged density in the pore, 
which in a bulk phase equals the usual reduced density, 
TO* = mj (ea^)^!'^ as the reduced dipole moment, and 

— ^w/f- for describing attractive walls (Eq. ([7])). 

For calculating the expansion of the dipole-dipole in- 
teraction Mayer function we used a second order approx- 
imation. This choice restricts the expansion coefficients 
aim{z) of the orientational distribution function to those 
with I < 2 and thus ItoI < 2. 



A. Phase diagrams 

As a first step, on the basis of the corresponding bulk 
theory we have determined the bulk phase diagrams of 
Stockmayer fluids for different dipole moments. For com- 
parison and as a reference the corresponding bulk phase 
diagram is depicted in all figures discussing the phase 
diagrams in slabs. We note that the bulk phase dia- 
grams presented here differ numerically from those in 
Ref because there a fourth-order expansion with re- 
spect to the corresponding Mayer functions has been used 
as well as a temperature dependent hard sphere diameter 
for the reference system. Extending our present analysis 
of spatially inhomogeneous system also to fourth order 
in that expansion is in principle possible but requires 
unreasonably large technical efforts. The bulk phase di- 
agrams within the present approximation exhibit only 
minor quantitative differences from those in RefJ^. 

For the dipole moment to* = 1.5 and the wall distance 
L/a = 10 Fig. 1 presents our numerical results for the 
phase diagrams in the chemical potential ~ temperature 
and the density - temperature planes. Figure 1 com- 
pares the bulk phase diagrams with those of the slab for 
repulsive or attractive walls. For describing the attrac- 
tive walls (Eq.([7])) here and in the following we choose 
= 1/2. Figure 1(a) shows that below the triple point 
temperature Ttp the isotropic vapor coexists with the 
ferromagnetic liquid. The confinement by repulsive walls 
shifts the chemical potential ~ temperature phase dia- 
grams towards higher chemical potentials relative to the 
bulk data. In the case of attractive walls the shift is 
smaller and in the opposite direction. This figure also 
shows that for the chosen parameters the Kelvin equa- 
tion provides a good estimate of these shifts. Between the 
temperatures Ttp and Tcp of the triple point tempera- 
ture and of the liquid- vapor critical point, respectively, 
there are three possible phases: the isotropic vapor, the 
isotropic liquid, and the ferromagnetic fluid. The first- 
order phase transition between the isotropic liquid and 
ferromagnetic fluid turns into a second-order phase tran- 
sition above the tricritical temperature Ttcp- (For large 
p* and at high p* orientationally disordered or ferro- 
magnetic solid phases appear—, which are not considered 
here.) We flnd that, compared with the bulk, the tem- 
perature range for liquid-vapor and liquid - ferromag- 
netic fluid coexistence is narrowed by the conflnement, 
both for repulsive and attractive walls. According to Fig. 
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FIG. 1: Phase diagrams of bulk and confined Stockmaycr flu- 
ids for m* = 1.5 and Lja = 10 in (a) the chemical potential - 
temperature and (b) in the averaged density - temperature plane; 
P* = -[J- d-^pi^) is the spatially averaged overall number den- 

sity. In (a) solid (bulk phase), dashed (attractive wall, = 1/2), 
and dash-dotted (repulsive wall) lines are lines of first-order phase 
transitions. In (b) these curves enclose two-phase regions. The 
dotted lines CL denote lines of second-order phase transitions. CP, 
TCP, and TP denote the liquid-vapor critical, the tricritical, and 
the triple point, respectively. Three-phase coexistence at TP is as- 
sociated with and visible by breaks in slope of the corresponding 
first-order phase boundaries. The symbols + and * indicate phase 
equilibrium points predicted by the corresponding ft'elvin equation 
for attractive and repulsive walls, respectively. 



1(b) at low temperatures (T < Ttp) the high density fer- 
romagnetic fluid coexists with the low density isotropic 
vapor. Above the triple point temperature Ttp, at low 
and medium densities two isotropic fluids coexist, be- 
coming identical above the liquid-vapor critical tempera- 
ture Tcp. At higher densities, the isotropic liquid (with 



FIG. 2: Phase diagrams of bulk and confined Stockmaycr fluids (as 
in Fig. 1) for rrt* = 1.5, L/a = 4, and ej, = 1/2 in (a) the chemical 
potential - temperature and (b) averaged density - temperature 
plane. 

the lower density) and the ferromagnetic fluid (with the 
higher density) are separated by first-order phase tran- 
sitions which turn into second-order phase transitions 
above the tricritical temperature Ttcp- The critical line 
CL of second-order phase transitions divides the liquid- 
like thermodynamic states into an isotropic liquid phase 
and into a ferromagnetic fluid phase. In comparison with 
the bulk coexistence curves the shift of the first-order co- 
existence region to lower densities is more significant for 
repulsive than for attractive walls. The tricritical liquid 
densities are lowered by the confinement in comparison 
with the corresponding bulk ones whereas the tricritical 
vapor densities are increased. The liquid-vapor critical 
density is lowered due to the confinement by repulsive 
walls and slightly increased in the case of attractive walls. 
At a given temperature both confinements promote the 
formation of the ferromagnetic phase at lower densities. 
These findings are in qualitative agreement with the re- 
sults of Gramzow and Klapp*^, but the shift of the coexis- 
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FIG. 3: Phase diagrams of bulk and confined Stockmayer fiuids 
for m* = 2, L/a = 10, and ej, = 1/2 in (a) the chemical potential 
— temperature and (b) averaged density - temperature plane. The 
line codes and the meaning of the symbols are the same as in Figs. 
1 and 2. For this dipole strength the liquid-vapor critical point 
has disappeared both in the bulk and in the confined system. In 
(a) the attractive wall confinement barely shifts the phase bound- 
ary; nonetheless this shifted boundary crosses the bulk one twice. 
Within the Kelvin equation for the shift those crossings correspond 
to zeroes (i) and (ii) of the independently calculated surface ten- 
sion difference "fwj'f — lw,if (see the inset) which correspond to a 
contact angle Qff= 90° (see the main text). 



tence curves of confined systems relative to the bulk ones 
is more significant than predicted by their Iromogcneous 
local density theory. As shown in Fig. 2, a stronger con- 
finement {L/a = 4) does not change the topology of the 
phase diagram. However, in this case the shift of the 
phase boundaries relative to the bulk ones are more pro- 
nounced. The temperature interval Tcp — Ttp, within 
which liquid- vapor coexistence is thermodynamically sta- 
ble, becomes narrower with increasing confinement. All 
critical temperatures (Tcp, Ttcp, Ttp) arc lowered upon 



increasing the confinement. (The lowering of the liquid- 
vapor critical temperature for Stockmayer fluids with in- 
creasing confinement is in a qualitative aggrement with 
the Gibbs ensemble Monte Carlo simulation results of 
Richardi et a&.) For attractive walls the critical den- 
sity of the liquid-vapor equilibria does not change signifi- 
cantly whereas it is reduced for repulsive walls. Compar- 
ing our results for m* = 1.5 and L/a ~ A with those of 
Gramzow and Klapp^ we note that, in contrast to their 
findings, in the chemical potential - temperature plane 
the first-order phase boundaries of the bulk and of the 
confined systems do not intersect. 

For L/a ~ 10, increasing the dipole moment from 
m* = 1.5 to TO* = 2 changes the corresponding phase 
diagrams shown in Fig. 1 considerably and leads to the 
phase diagram displayed in Fig 3. This increase of the 
dipole moment causes the disappearance of the liquid- 
vapor transitions for both the bulk and the confined sys- 
tems. According to Fig. 3(a) for this dipole moment the 
phase boundaries of the bulk and of the confined system 
with attractive walls do intersect. On the basis of the 
Kelvin equation (Eq. ((39)) ) this multiple crossing of the 
^rromagnetic .^uid - isotropic ^uid (/ = //, // = if) 
phase boundaries for the bulk and for the confined system 
implies a non-monotonic temperature dependence of the 
surface tension difference Jwjf — Iw.if and vice versa; 
according to the full line in Fig. 3(b) Pff ~ Pif > 
varies monotonically as function of temperature. (For 
smaller dipole moments a similar crossing behavior has 
been found by Gramzow and Klapp^.) The inset in Fig. 
3(a) shows the independently calculated (Eq. (^0)) ) sur- 
face tension difference and confirms that the aforemen- 
tioned crossings of the bulk and film phase boundaries 
do coincide with the changes of sign of "fw.ff — IwSf at 
T(i) and T(^ny Due to the second equation in Eq. 
the Kelvin equation turns into 



A^^^ 2 7//,^/ cos 6// 
L Pff - Pzf 



(41) 



Since jffjf is expected to decrease monotonically upon 
increasing temperature, these crossings amount to a non- 
monotonic temperature dependence of the contact angle 
of the ferromagnetic fluid around 90°. (For 8// = 90° 
the leading behavior of Ap is given by the contribution 
^ 1/L\) 

In the case of confinement with repulsive walls the fer- 
romagnetic fluid - isotropic fluid phase boundary is gen- 
erally shifted to higher values of the chemical potential as 
compared with its bulk counterpart. Confinement lowers 
the tricritical temperatures Ttcp and there is no sig- 
nificant difference between the tricritical temperatures 
for systems with attractive and repulsive walls. As for 
TO* = 1.5 (Fig. 1(a)), also for to* = 2 and L/a = 10 the 
estimates obtained from the Kelvin equation are in good 
agreement with the full phase equilibrium calculations. 
Figure 3(b) shows how the two-phase region shrinks upon 
confinement. If the liquid- vapor critical point has disap- 
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FIG. 4: Phase diagrams of bulk and confined Stockmayer fiuids 
for m* = 2, L/ct = 4, and ej, = 1/2 in (a) the chemical potential 
— temperature and (b) averaged density - temperature plane. The 
line codes and the meaning of the symbols are the same as in Figs. 
1-3. In (a) the crossings of the shifted phase boundary for the 
attractive wall confinement with the bulk phase boundary is more 
clearly visible than for Lja = 10 shown in Fig. 3(a). 

peared in the bulk phase diagram, this holds also for 
the confined system, independent of the character of the 
wall. For m* = 2, upon decreasing the wall separation 
from L/a = 10 to L/a = A the topology of the phase di- 
agrams does not change (compare Figs. 3 and 4). Figure 
4(a) shows that for repulsive walls the stronger confine- 
ment leads to a larger shift of the phase boundary to 
higher values of the chemical potentials. For attractive 
walls this shift remains much smaller even for this nar- 
row slab with L/cr ~ A. Different from the case m* ~ 1.5 
(see Fig. 2(a)), for m* = 2 two crossings of the bulk 
coexistence curve with the phase boundary for attractive 
walls are clearly visible. Their interpretation in terms 
of the Kelvin equation is less valid than for Fig. 3(a) 
because due to the small value L/a = A subdominant 
terms in the expression for A/i become relevant. Fig- 



ure 4(b) shows that enhancing the confinement causes 
a further shrinking of the two-phase coexistence region. 
The tricritical temperatures Ttcp and densities prcp 
decrease upon reducing the wall separation. This means 
that for smaller distances between the walls the ferro- 
magnetic fluid - isotropic fluid phase coexistence occurs 
within a narrower range of the thermodynamic variables. 

In a recent publication Trasca and Klappi have stud- 
ied inter alia the second-order phase transition of strongly 
coupled dipolar fluids conflned to narrow slit pores. They 
performed Monte Carlo simulations using purely repul- 
sive wall confinements and the dipolar soft sphere model 
for the particles. In confined systems, they have found 
that, at a given temperature, with the decrease of the wall 
separation the paramagnetic-ferromagnetic phase transi- 
tion density (averaged across the slit pore) decreases rel- 
ative to the corresponding bulk one. They found that 
the direction of this shift is inconsistent with their very 
simple mean-field theoretical predictions (sec Eq. (3a) in 
Refji). Analyzing the shifts of the second-order transi- 
tion lines for the repulsive wall confinements (see Figs. 
3(b) and 4(b)) we find that within the framework of our 
DFT at a given temperature the critical averaged density 
decreases upon the decrease of the wall separation, which 
is in qualitative aggrement with the simulation findings of 
Trasca and Klapp*^. (There is no possibility for a quanti- 
tative comparison because these Monte Carlo simulations 
were carried out for a reduced dipole moment m* ~ 3; 
for such large values of m* the quantitative reliability of 
the present DFT would be reduced anyhow.) 



B. Structural properties 

In the following the structural properties of the coex- 
isting ferromagnetic fluid and isotropic gas phases are 
discussed. Figure 5 shows the corresponding results for 
reduced dipole moment m* = 1.5 and wall separation 
L/a = 10 at the reduced temperature T* = 1.15. At this 
low temperature the ferromagnetic fluid can be denoted 
as a liquid (see Fig. 1(b)). (As a hint we note that for 
a given temperature the chemical potential at two-phase 
coexistence is different for attractive and repulsive walls, 
see Fig 1(a).) Figure 5(a) shows that the confined fer- 
romagnetic liquid phase is strongly structured both for 
attractive and repulsive walls. In the case of confine- 
ment by attractive walls the contact value ~ 19.8 of 
the reduced density (i.e., at \z\ ~ (i — cr)/2) is very high 
(not displayed on the scale of the figure) which refers 
to a strong adsorption on the walls. The density pro- 
files of the coexisting isotropic gases are displayed in Fig. 
5(b). For the attractive wall case the contact value of 
the density is ~ 0.55 which indicates a visible spatial 
inhomogeneity in the vicinity of the walls. However, in 
the case of repulsive walls, at this low temperature the 
density of the vapor phase does not show any detectable 
inhomogenity. The behavior of the preferential orienta- 
tion a2o{z) shown in Fig. 5(c) tells that in the liquid 
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FIG. 5: Profiles of the number density p*{z) ((a) and (b)), pref- 
erential orientation 020(2) of tlie dipole axes (Eq. II14I I) ((c) and 
(d)), ferromagnetic order parameter axy(z) (Eq. I I33I I) (c), and 
in-plane magnetization M*y{z) = m* p* (z)axy(z) (f) of coexist- 
ing ferromagnetic liquid and isotropic vapor phases for m* = 1.5, 
L/cr = 10, and T* = 1.15, confined by repulsive ( — ) or attractive 
{■■■■,e^ = 1/2) walls. 

phase the dipole axes are preferentially oriented parallel 
to the walls throughout the pore for both types of con- 
finements. The ferromagnetic ordering of the particles in 
the liquid phase is confirmed by the variation of the order 
parameter axy{z) (Eq. ((33|) ) which can be seen in Fig. 
5(e) which describes the net orientation of the dipoles 
parallel to the walls, keeping in mind that for the orien- 
tational order parameter we find aio{z) = {Mz ~ 0, 
see Eq. ([5^ ) throughout of the pore. We note, that 
for all confined fluid phases studied here we found that 
aio(z) = 0, i.e., throughout the pores the z component 
of the magnetization is zero. Figure 5(f) displays the 
corresponding magnetization Mxy{z) (Eq. ([Ml) ), which 
shows strong spatial inhomogeneities for both types of 
confinements induced by the structure of p{z). (For the 
attractive walls the contact value of the magnetization is 
Mxy^w — 25.7.) For the vapor phase we find axy{z) = 
for both attractive and repulsive confinements, which to- 
gether with Mz = implies that the vapor phase is an 
isotropic phase. However, in Fig. 5(d) the order param- 
eter a2o{z) shows that close to the walls there is a slight 
ordering in the vapor phases as well. There the order- 
ing is more pronounced for the confinement by attrac- 



FIG. 6: Same as Fig. 5 for L/ct = 4. 



five walls than by repulsive walls. In the vapor phase, 
close to the attractive walls the dipole axes arc prefer- 
entially ordered parallel to the walls (a2o < 0), but the 
amplitude of the orientation is much smaller than in the 
coexisting liquid phase (Fig. 5(d)). Between this layer 
of parallel orientation and the orientationally disordered 
interior of the pore (i.e., a2o{z) is close to zero), there is 
a layer with preferential orientation perpendicular to the 
wall (q;2o > 0, Fig. 5(d)). With respect to the orienta- 
tional order, our DFT results are consistent with those 
of Gramzow and Klapp^ but they show stronger spatial 
inhomogenities which is borne out by using a more so- 
phisticated free energy functional for the hard sphere ref- 
erence system. For m* = 1.5 and T* ~ 1.15, upon de- 
creasing the wall separation from L/a = 10 to L/cr = 4 
the density and orientational distributions change signif- 
icantly. Figure 6(a) shows that for both attractive and 
repulsive confinements the ferromagnetic liquid phases 
become more structured. Moreover, the confinement by 
attractive walls induces more ordered structures than by 
the corresponding repulsive ones. The density profiles for 
the coexisting isotropic vapor phase are displayed in Fig. 
6(b), telling that only the attractive walls induce a spatial 
inhomogeneity in the density distribution of the isotropic 
vapor. Figure 6(c) shows that in the liquid phase the ex- 
tent of orientation of the dipole axes parallel to the walls 
is comparable with that in the wider {L/a = 10) pore 
but it varies more strongly across the pore. According 
to Fig. 6(d) the same conclusion holds for the structural 
properties of the confined vapor phase; since the vapor 
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FIG. 7: Same as Fig. 5 for increased temperature (T* = 1.6) and 
stronger dipole moment (m* = 2). 



FIG. 8: Same as Fig. 7 for L/a = 4, 



phase is isotropic this shows that the preferential order- 
ing parallel and perpendicular to the walls does not hinge 
on the formation of ferromagnetic order. From Fig. 6(e) 
one infers that the in-plane ferromagnetic order param- 
eter of the liquid phase for L/a = A is comparable with 
that for L/a = 10. However, the thinner pore exhibits 
a stronger magnetization and more pronounced spatial 
variations (Fig. 6(f)). 

Figure 7 displays the structural properties for the wide 
pore L/a = 10 at the higher temperature T* — 1.6 and 
for an increased dipole moment m* = 2. Compared with 
Fig. 5 (T* = 1.15, m* = 1.5) there are no qualitative 
differences. The quantitative differences concern the va- 
por density, which is lower due to the higher temperature 
(Fig. 7(b)), and the magnetization of the liquid phase, 
which is higher due to the increased dipole strength (Fig. 
7(f)). The occurrence of preferential ordering parallel to 
the walls in the vapor phase (0120 (z) < 0, Fig. 7(d)) is in 
line with the findings of Gramzow and Klapp^; however 
their approach does not render the preferential ordering 
perpendicular to the walls {a2o{z) > 0) in the subsequent 
layer towards the interior of the pore. But this prefer- 
ential orientation is about a factor of ten smaller (Fig. 
7(d)) as compared with the case shown in Fig. 5(d). 

For TO* = 2 and T* = 1.6, decreasing the wall sepa- 
ration from L/a = 10 to L/a = 4 leads to a significant 
increase of both the liquid and the vapor densities (Figs. 
8(a) and (b)) as well as to more pronounced density oscil- 
lation. The preferential oricntional ordering is similar for 



the narrow (Figs. 7(c)-(e)) and the wide (Figs. 8(c)-(e)) 
pore. The magnetization is stronger in the more confined 
pore (Fig. 8(f)); this and the more pronounced spatial 
variation is caused by the corresponding behavior of the 
underlying number density (Fig. 8(a)). The vapor phase 
does not acquire a ferromagnetic order even in the narrow 
pore, i.e., axy{z) = for the vapor phase. 



VIII. SUMMARY 

In the present study of the phase behavior and struc- 
tural ordering of confined Stockmayer fluids (Eq. ^) 
the following main results have been obtained. 

1. We have applied an extension of the modified mean- 
field density functional theory to determine the phase 
behavior and the structure of Stockmayer fluids in slit- 
like pores formed by either attractive (Eq. ([7])) or re- 
pulsive (Eq. ^) walls at wall separations L/a ~ 10 
and L/a ^ i, for reduced dipole moments to* = 1.5 and 
TO* = 2, and for a reduced wall strength = 1/2. This 
system exhibits three distinct fluid phases: isotropic va- 
por, isotropic liquid, and ferromagnetic fluid. 

2. For the reduced dipole moment to* = 1.5 wc have 
found isotropic liquid - isotropic vapor, ferromagnetic 
fluid - isotropic vapor, and ferromagnetic fluid - isotropic 
liquid first-order phase separations. The confinement by 
repulsive walls shifts the phase boundaries in the chem- 
ical potential - temperature plane towards higher chem- 
ical potentials while attractive walls lead to a shift to- 



12 



wards lower chemical potentials relative to the corre- 
sponding bulk data. The decrease of the wall separation 
from i/(T = 10 to L/cr = 4 docs not change the topology 
of the phase diagrams, but in this case the shift of the 
phase boundaries relative to the bulk ones are more pro- 
nounced (see Figs. 1 and 2). The increase of the dipolc 
moment from m* = 1.5 to m* = 2 erases the isotropic 
liquid - isotropic vapor phase transitions for both the 
confined and the bulk systems. We have found that for 
this stronger dipole moment only the phase boundary of 
the confined system with repulsive walls is significantly 
shifted relative to the bulk one (see Figs. 2 and 3). The 
phase boundaries of the bulk and of the confined system 
with attractive walls do intersect without a significant 
shift of the phase boundaries of the confined system rel- 
ative to the bulk one. 

3. We have shown that in the temperature - chemical 
potential plane the Kelvin equation provides a good es- 
timate of the confinement induced shift of the chemical 
potential for isotropic vapor - ferromagnetic fiuid (Fig. 
1(a)) and isotropic fluid - ferromagnetic fluid (Fig. 3(a)) 
phase coexistence. The multiple crossings of the isotropic 
fluid - ferromagnetic fluid phase boundaries in the tem- 
perature - chemical potential plane correspond to a non- 
monotonic temperature dependence of the contact angle 
of a drop of ferromagnetic fluid on a single wall. 

4. On the basis of our DFT approach we have found that 
the number density p{z) of the confined liquidlike phases 
(coexisting with the corresponding vaporlike phases) for 
both dipole moments (m* = 1.5 and m* = 2) and for 
both confinements {L/a = 10 and L/a = 4), with ei- 
ther attractive or repulsive walls, are strongly structured 
(see Figs. 5(a)-8(a)). The figures for the order parameter 
a2o{z) (Eq. (HU) for preferential orientation of the dipole 
axes show that in the liquidlike phases the dipole axes 
are preferentially oriented parallel to the walls through- 
out the pore for both attractive and repulsive walls as 
well as for both sizes of confinements (see Figs. 5(c)- 
8(c)). The ferromagnetic ordering of the dipoles in the 
confined liquidlike phases is displayed by the variation 
of the ferromagnetic order parameter axy (see Eq. p3p 
and Figs. 5(e)-8(e)) which, together with the result 
Mz{z) = 0, describes the net orientation of the dipoles 
parallel to the walls. The corresponding magnetizations 
Mxy{z) ~ mp[z)axy{z) are displayed in Figs. 5(f)-8(f). 
Since axy{z) = and q;io(-z) = throughout the pore for 
the vaporlike phases (at coexistence with the ferromag- 
netic liquidlike phases), they are completely isotropic. 
Attractive wall confinements give rise to strong adsorp- 
tion at the walls even in the vapor like phases (see Figs. 
5(b)-8(b)). In the vapor like phases, close to the attrac- 
tive walls a slight preferential parallel ordering has been 
found. Between the layer of this preferential parallel ori- 
entation and the orientationally disordered interior of the 
pore a layer with preferential perpendicular orientation 
has been detected (see Figs. 5(d)-8(d)). 
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APPENDIX A: FUNDAMENTAL-MEASURE 
THEORY 

One of the main limitations of the original version of 
fundamental-measure theory is that the underlying bulk 
fluid equation of state reduces to the Percus-Yevick com- 
pressibility equation. As it is well known, within this 
approximation the contact value of the density profile of 
a one-component hard-sphere fluid at a planar wall is sig- 
nificantly overestimated at high bulk densities. In order 
to improve the quantitative accuracy and to maintain in- 
ternal consistency, following Roth et ali^i and Tarazona— 
here we apply a fundamental-measure theory based on 
the Carnahan-Starling equation of state. Within this ap- 
proach the hard-sphere excess free energy density func- 
tional is given in terms of weighted densities as 



fhs 



<f({"a(^)}), 



(A42) 



where 



$ = -noln(l - na) 



ni7i2 - ni-n2 

1 - 713 



,/3 Q , ria + (1 - "3)^ln(l - ng) 

-f - 3n2n2-n2 — — j- . A43 

3677713(1 — n^y 

In Refi^ it is shown that the weighted densities in the 
present slab geometry are given as 

na{z)= dz'p{z + z')w^"\z'), (A44) 

J -a/2 

where the reduced weight functions as functions of 
z are 

2 



a 



^(1) = ^e^, w(2) = 27rze^. (A45) 

The chemical potential functional (first-order direct cor- 
relation function) of the reference hard sphere system is 



5Fhsmz)}] 
5p{z) 



t r\z''-^^.i-\z' ^ z), (A46) 
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where the corresponding derivatives are 



where 



a* 

dno 



$1 = 

<I>2 = 



712 



dni 1 — 713 ' 
9$ ni 



+ {nl - n2n2) 



dn2 1 — 713 
(713 + (1- 713)2 hi(l-n3)) 



127rn2(l - 113)2 



$3 



Ho 



(9713 1-713 

ni7i2 - nin2 , 



(lij - 37i2n2n2) 



(1-^3)2 

/ 2 + 713(713 -5) ln(l - 713) 
Iv367r7i2(l- 113)3 187r7i3 
9$ n2 



(f>4 



(9ni 1 — 713 ' 
9$ ni 



8X12 1-713 



713 + (1- 713)2 ln(l- 713) 

""2n2 T 27i T2 ■ (^47 

D7r7i3(l — 713)"^ 



In the summation (Eq. (|A46[) ) a = 4, 5 referes to the 
vectorial weight hmctions w'^^ and w^^^ respectively. 



X= cos 6^12 



The coefficients ^Ijl^^j^ma symmetric, 



Vi?2 + Z2 



(B49) 



4:in^.™.(^)=A^:Uh™,w. (B50) 

so that the functions 

are symmetric, too, which simplifies Eqs. (|36p and (j37p . 
We note that the integral in Eq. (pS)) splits into two parts 
due to the Hcavisidc function 0(i?2 ^2 — 0-2): 

/"OO 

Jo 

= e{\z\-a) / dRg{R,z) 
Jo 



+Q{<J-\z\) / ^ di?5(i?,z) 



(B52) 



APPENDIX B: CALCULATION OF THE 



COEFFICIENTS 



h^2 



Solving the set of integral equations ([55)1 and ((57|) (see 
also Eqs. (|C53p and (jC54p ) requires to determine the 

(z) (see Eq. ^7} ) which are based on 



functions 



Iirnil2m2 



the coefficients A 



Iiniil2ni2 



(x) (see Eqs. and 
By using Eq. we have calculated these coefficients 

and in the following we provide all coefficients which are 
nonzero for i, li, I2 = 0, 1, 2: 



_ .(1) _ Stt^ 
^0000 ~ ' ^1010 ~ 2 


-(3x2 _ 


-1), 
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.(2) _ 167r2 (2) _ SVStt^ 
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^2020 ~ g W-^ 
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1(2) 
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-2121 - 15 
l(2) _ 



(-18x^ + 15x2 
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2222 



6x^ 



1), 

- 1), (B48) 



APPENDIX C: CALCULATION OF THE 
DENSITY AND ORIENTATIONAL PROFILES 



Equations (|36p and ([57)1 lead to coupled integral equa- 
tions for p{z) and a(z,a;) (or q;;„i(2), see Eq. (fT3|) : here 
we suppress their dependences on T and fi): 
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p(z) = exp ~ PfJ-h^s^iipiz)}] — / duja{z,uj)ln[iTTa{z,Lu)] 



l'L/2 

-2/3 ^ ahmi(^) X! / dz'p{z')ui^,n^i^,n.^{z- z')ai.^^.^{z')+ Py.- (iu.^{z) 
Elementary calculations lead from Eq. ([57]) to the following equation: 



(C53) 



a(z, uj) 



L 11 

/ exp [-2/3 /2 Eh^mi (^') T.l^,r,^^ "iimi/2™2 - ('^)] 



(C54) 



which fulfills the normalization requirement in Eq. ([5]). Using the expansion of a(z, w) in terms of spherical harmonics 
(Eq. (fT5|) ) one obtains for the coefficients a/„i(z) (Eq. P^ ): 



, , /f^'^^;tmi(^)exp[-2/3/„^,2'^zV(z')E;2,rn,"i2m2(^')E;3,rn3"'2™2i,3m3(2-2')y;;„3(t^)] , 

ahmi(z) = • (055) 

/ exp [-2/? /_^/2 dzV(z') Ei,,m2 "'2^2 (^') Ei3,m3 Mi2m2i3m3 " ^O^iarng (^)] 



Equations ((C53)) and ((C55)) are a set of coupled integral entropy term must be expressed in terms of the coefR- 

equations which are solved iteratively for the density pro- cients aim(-z)- Using Eq. (jC54[) . after some calculations 

file p(z) and the coefficients aimiz) of the orientational one arrives at 
distribution. To this end. in Eq. (jC53p the orientational 



c?a;a(z, w) ln[47rQ;(z, w)] = 
ln[47r/C(z)] - 2/3 ^ ai^m^z) \ dz p{z')ui^rn.^i^m^{z - z')ai^Tn^{z'), 



/i ,mi /2 ,^2 



where 



r rL/2 

C(z) = / dojcxp -2/3 / dz'p{z') V a,,,„,,(z') V u/.^ii^m^ (2; - z')y;*„,^ (w) 



^2 -''7^2 



(C56) 



(C57) 
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